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ABSTRACT 



Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful 
probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the 
time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday 
rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total 
electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of 
this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday 
rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation 
measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations — either intrinsic or due to 
the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the 
ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, 
as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, i 
such as the SKA and its pathfinders. 



Key words. Polarization - Techniques: polarimetric 



1. Introduction 

In recent years, low-frequency (v Q b s < 300 MHz) radio astron- 
omy has undergone a revival due to the construction of the Giant 
Metrewave Radio Telescope (Swamp 1991 GMRT;) and mod- 
ern aperture array radio telescopes such as the LOw Frequency 
ARray (LOFAR; Stappers et al. 2011; van Haarlem et al 2013, 
submitted), the Long- Wavelength Array (Kassim et al. 2010 



LWA;), and the Murchison Widefield Array (Mitchell et al 



2010| MWA;). The first phase of the Square Kilometre Array 
^ Garrett et al.||2010| SKA;) is also planned to feature a large 
number of low-frequency antennas, operating at ~ 70 - 450 
MHz. These telescopes open new scientific horizons in the area 
of low-frequency radio astronomy, including the determination 
of high-precision Faraday rotation measures (RMs). 

Faraday rotation causes the intrinsic polarization angle (xo) 
of a signal to rotate as it propagates through a magneto-ionic 
mediunr] The observed polarization angle of a point source can 
be defined a^] 



X = X0 + (<Pion + 0ISM + 0IGM)/1 



(1) 



where x denotes the observed polarization angle in radians and 
A the observing wavelength in meters. <p is the Faraday depth 
in rad irT 2 of a given intervening magneto-ionic medium. In 
Eq. [T] three intervening media produce Faraday rotation: the 
ionosphere (ion), the Galactic interstellar medium (ISM), and 
the inter-galactic medium (IGM; assuming the source is extra- 
Galactic). Low-frequency observations are particularly affected 
because Faraday rotation increases quadratically with wave- 
length. 

For the case of a single polarized source positioned behind 
one or more magneto-ionic media that are not emitting polarised 
radiation, the Faraday depth of the source is equivalent to its 
RM (Sok oloff et al.|1998 e.g.,). Nonetheless, here we will use 
the more generic term Faraday depth. Following Brentjens & de 
Bruyn (2005) we define Faraday depth as: 



0(1) = 0.81 
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»^ SO 



■observer 



h p B ■ dlradm 



(2) 



where n e and B are the electron density (cm -3 ) and mag- 
netic field (pG) integrated along the line-of-sight (LOS) to the 
source and dl is the infinitesimal path length in pc. A mag- 
netic field pointing towards/away from the observer gives a pos- 
itive/negative Faraday depth. 

The electron density in the ionosphere dictates the lowest 
frequency observable from the ground, approximately 10 MHz. 
Above this frequency, the ionosphere affects signals in three 
main ways: i) differential phase delays, ii) Faraday rotation, and 
iii) absorption in the High Frequency band (HF; 3-30 MHz) and 
the low-end of the Very High Frequency band (VHF; 30-300 
MHz) due to the presence of the so-called "D-region" in the 
daytime. Assuming a typical observing frequency of 150 MHz 
(LOFAR high band) and an ionospheric Faraday depth of 1 rad 
m -2 , the additional rotation of the polarization angle imparted 
by the ionosphere will be ~ 228.9°. Although the rotation of the 
polarization angle is less pronounced at higher frequencies, the 
Faraday depth of the source will still be systematically affected 
by the ionosphere. Due to the direction of the geomagnetic field, 



1 A magneto-ionic medium is made up of ionized gas and magnetic 
fields (e.g., the ionosphere). 

2 First-order approximation for high frequencies. In the low- 
frequency band additional terms, not shown, may also become signifi- 
cant. 



ionospheric Faraday rotation has a positive or negative contribu- 
tion to the total Faraday depth of a source when observing from 
the northern or southern hemispheres, respectively. For instance, 
the contribution from the ionosphere should be corrected for in 
order to derive reliable Faraday depths due to the ISM alone 
when observing Galactic pulsars. This is particularly important 
for pulsars that are relatively nearby and/or located above the 
Galactic plane because the magnitude of the ionospheric Faraday 
depth can be a significant fraction of, or even greater than, the 
total observed Faraday depth. 

Calibrating for ionospheric Faraday rotation is complicated 
because the free electron content of the ionosphere varies de- 
pending on the time of day, season, level of solar activity, 
and LOS. The ionosphere changes on timescales that are often 
shorter than the length of an observation; e.g., Brentjens] (2008) 
andlPizzo et al. (2011 1 report Faraday depth variations in polar- 
ized point sources of a few rad m -2 in 12 hours. Calibrating for 
the ionosphere is therefore critical for comparing Faraday depths 
of the same source at multiple epochs. The time-dependence of 
ionospheric Faraday rotation can wash out the linear polariza- 
tion when averaging over multi-hour time intervals at long wave- 
lengths. 

Faraday depth measurements can be used to map the struc- 



ture of the Galactic magnetic field (GMF) using pulsars (Han 
|et al.|200rj Noutsos et al.|2008] e.g., ) and extragalactic sources 
( |Brown et al.||2007[ |Van Eck et al.||2011) . Knowledge of the 
magnitude and structure of the GMF is key for understand- 
ing deflection of high-energy cosmic rays, star-forming regions, 
instability-generated turbulence, pressure on ionized gas and 
the transport of heat, angular momentum and energy from cos- 
mic rays. Monitoring Faraday rotation variations over time also 
yields insights into the polarization modes of pulsar emission 
and ISM magnetic field variations (Weisbe rg et al.|2004| e.g.,). 
The magnetic fields of other galaxies have also been the sub- 
ject of thorough investigation, although weak magnetic fields 
best detected at low frequencies remain relatively unexplored 
( Beck|2009 i. Detecting weak, coherent magnetic fields will pro- 
vide insight into how distant cosmic-rays originating in the discs 
of galaxies propagate within the halos and possibly into inter- 
galactic space. Furthermore, a deeper understanding of mag- 
netic fields in galaxy cluster halos and relics can also be gained. 
Planetary observations can be used to determine ionospheric 
Faraday depths, e.g., by observing and studying the nature of 
Jupiter bursts at ~ 20 MHz ( |Nigl et al.||2007[ ). Lastly, Faraday 
depth measurements may even provide a test for the nature of 
the accretion flow of the supermassive black hole at the Galactic 
Centre ( |Pang et al.|2011| e.g.,). 

Here we present ionFR, a code that models the ionospheric 
Faraday depth using publicly available global TEC maps and 
the latest geomagnetic field model. We demonstrate the ro- 
bustness of ionFR by comparing its modeled Faraday depths 
with the measured Faraday depths of pulsars from four LOFAR 
observing campaigns, as well as an observation of the pulsar 
PSR B 1937+21 with the Westerbork Synthesis Radio Telescope 
(WSRT). §2 describes previous studies and the methodology 
of the ionFR code. §3 presents modeled ionospheric Faraday 
depths that the software produces for varying levels of solar 
activity, including for the proposed SKA sites. §4 briefly in- 
troduces RM-synthesis, which was used to determine Faraday 
depths here. Comparisons between the ionFR-modeled iono- 
spheric Faraday depths and the observations are shown in §5. §6 
discusses these results in more detail. Conclusions are presented 
in §7. 
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2. Modeling the ionosphere with the ionFR code 

Here we describe the theoretical background and methodology 
of ionFR, which is written almost entirely in Python and is freely 
available to the community via sourcef orgqj The code returns 
a table containing the ionospheric TECj^] magnetic field magni- 
tude and RM along the requested LOS (see Fig. [4] for two ex- 
amples). The required input arguments are: right ascension and 
declination of the source, geographic coordinates of the observ- 
ing site, the date of observation and the type of TEC map to be 
used. 



2.1. Previous implementations 

There have been several ionospheric Faraday rotation models 
previously presented. Erickson et al. (2001) constructed a sim- 



pie model based on Global Positioning System (GPS) data. 
However, this model requires local GPS data, i.e., dual frequency 
GPS receivers installed at the telescope site(s). Afrai movich| 



et al. 



mode 



(2008) used a well-established empirical ionospheric 
, the International Reference Ionosphere (IRI). The IRI 
provides ionospheric parameters derived mostly from ground- 
based instruments (e.g., ionosondes) and some space-based in- 
struments (Bilitza & Reinisch 2008 ). The model was compared 
solely with GPS data and no comparison with radio astronom- 
ical data was presented. The ionFR software presented here 
is somewhat similar to the TEC0R task from the Astronomical 
Imaging Processing System ( |Greisen||2003| AIPS;). However, 
while TECOR requires exporting interferometric data into AIPS, 
ionFR is a standalone package. 



2.2. Ionospheric piercing point 

To calculate the ionospheric Faraday depth along the LOS, we 
assume a thin spherical shell surrounding the Earth (Fig. [T]>. The 
ionospheric Faraday depth is then calculated at the ionospheric 
piercing point (IPP). Under these assumptions and from Eq. [2] 
the ionospheric Faraday depth is defined as: 



2.6 x 10" 17 TEC L os B l0S rad m 2 , 



(3) 



where TEClos is the total electron content at the geographic co- 
ordinates of the IPP. TEClos is defined as: 



TEClos 



= fn t 



dim" 



(4) 



Blos is the geomagnetic field intensity in gauss at the IPP. To 
facilitate the estimation of the IPPs, the code assumes that the 
Earth is a sphere of radius R = 6371 km. 

Given the triangle defined by the points A, O, and IPP 
(Fig. [Ti, the value of the zenith angle (£') at the IPP is derived 
using the law of sines: 



sin(r) 



R 



R + h 



sin(<D, 



(5) 



where h is the altitude of the ionospheric thin shell, as specified 
in the ionospheric electron density files described in §2.3 In a 



similar fashion, the other geographic and topographic parame- 
ters at the IPP can be calculated. Spherical trigonometry is used 
to calculate the latitude, longitude, and azimuth (ip') of the IPP. 

3 http : //sourcef orge . net/pro j ects/ ionf arrot/ 

4 Typically measured in TEC units (TECU) 
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Fig. 1. Schematic representation of an astronomical signal pierc- 
ing the ionosphere. The solid blue curves delineate the iono- 
sphere. The thick red curve is the ionospheric thin shell approx- 
imation, and the thin black curve delineates the surface of the 
Earth. 

2.3. Ionospheric electron density 

Measurements of the ionospheric free electron content come 
from several sources. For example, the Center for Orbit 
Determination in Europe (CODE) offers global ionosphere maps 
(GIMs) in IONosphere map EXchange format (IONEX), avail- 
able via anonymous £trj^| IONEX files from CODE are derived 
from ~ 200 GPS sites of the International Global Navigation 
Satellite System Service (IGS) and other institutions. Figure [2] 
illustrates twelve maps obtained from the CODE IONEX file for 
April 11th, 2011. 

The IONEX files provide vertical total electron content 
(VTEC) values in a geographic grid (Ai on = 5°, Ai at . = 2.5°). 
VTEC is defined as the integral of free electrons in the iono- 
sphere along the zenith (£ = 0°) direction. The time resolution 
of the IONEX files provided by CODE is 2 hours. To increase 
the time resolution, ionFR creates new GIMs for every hour by 
using an interpolation scheme that takes the rotation of the Earth 
into consideration (Sch aer et al.|19 98, Eq. 3). The software in- 
terpolates the positional measurement grid to estimate the VTEC 
values at a given IPP for each hourly GIM. The interpolation 
uses a 4-point formula as in Fig. 1 of Schaer et al.| ( 1998). 

The final step in calculating the line-of-sight TEC is convert- 
ing the VTEC to the slant TEC (TEClos) as follows: 



TEClos 



VTEC 
cos(^')' 



(6) 



CODE also provides root-mean-square (RMS) VTEC maps that 
are geographically gridded in the same way as the VTEC maps. 
The uncertainties are calculated from these maps using Eq. [6] 
The 1-cr uncertainties in the RMS VTEC maps are typically be- 
tween 2-5 TECU. 

ionFR will be regularly updated to allow a greater selection 
of TEC map sources. Maps with higher spatial and temporal res- 
olution are desirable to trace ionospheric variations on shorter 



ftp : //ftp . unibe . ch/aiub/CODE/ 
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Fig. 2. GIMs representing the VTEC across the globe for April 1 1th, 201 1 (the date of the first LOFAR observing campaign, see §5) 
obtained courtesy of CODE. The maps range from minimum (blue) to maximum (red) VTEC values of 0.0 - 87.2 TECU (1 TECU 
= 10 16 electrons/m 2 ). The triangles indicate the location of the LOFAR core stations in the Netherlands, the squares mark the SKA 
core sites in South Africa and Western Australia, and the circles indicate the site of the GMRT. 



timescales and smaller spatial scales. For European telescopes, 
ionFR can also use TEC maps from the Royal Observatory of 
Belgium (ROB0), which are derived from GPS data from a per- 
manent European network. These TEC maps are more finely 
gridded than those from CODE (A lon . = 0.5°, A lat . = 0.5°) and 
have 15-minute time resolution (Fig J3]l. They are also now pub- 
licly available via anonymous ftrFTand are being produced 
since the beginning of 2012. Comparisons of ionFR-modeled 
Faraday depth based on the CODE/ROB maps are discussed in 
§5.1 (see Figs. |]and|9]). 

2.4. The geomagnetic field 

The Earth's magnetic field is calculated using the eleventh 
generation of the International Geomagnetic Reference Field 
(IGRF11; |Finlay et aL] |20T0] > released in December 2009. 
The IGRF is derived by the International Association of 
Geomagnetism and Aeronomy (IAGA) every five years and is 
available to downloacj^] The IGRF1 1 is described as the negative 
gradient of a scalar potential, B= - VV, which is a finite series of 
spherical harmonics. ionFR calls the IGRF1 1 to deliver the vec- 
tor components of the geomagnetic field at the IPR These point 



| http://gnss.be/Atm ospheric_Maps/ ionospheri c_maps . php| 
[ftp : //gnss . oma . be/gnss/products/I0MEX7| 
http : / /www . ngdc . noaa . gov/IAGA/vmod/igr f . html 





Leys m.^E 



Fig. 3. The VTEC across Europe for March 23rd, 2012 (the date 
of the third LOFAR campaign, see §5) at 00:00 UT, obtained 
courtesy of ROB. The square indicates the LOFAR core stations 
and the triangles represent the locations of the international sta- 
tions. 



towards the north (X), east (Y), and radially towards the center 
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of the Earth (Z). The XYZ coordinates are local (orthogonal) co- 
ordinates, i.e. X is not pointing in the direction of global north 
but to the local (tangential) north on the sky. The total magnetic 
field along the LOS at the IPP is then estimated as follows: 



B L0S = Zcos(£') + Fsin(Osin(</>') + Xsin(£')cos(<p'). (7) 



2.5. Error propagation 

The fractional uncertainties on Blos compared with their central 
values (ct Blos /Blos) are much smaller than the fractional uncer- 
tainties in the slant TEC (ct TE c los /TEClos)- Consequently, to 
determine the uncertainties on 0; on , only the RMS slant TEC val- 
ues are used. This results in uncertainties of 0.1 - 0.3 rad irT 2 in 
</>ion using the CODE global TEC maps. Assuming RMS values 
of 0.5 TECU for the ROB European TEC maps results in smaller 
uncertainties of 0.03 - 0.06 rad irT 2 . 



3. Ionospheric RM variation 

Here we present some examples of ionFR-modeled ionospheric 
Faraday depths using TEC data from CODE. We illustrate both 
the 0; on variation within a day (for two epochs with differing 
levels of solar activity) as well as the longer-term variations with 
season and solar activity. 

The modeled ionospheric Faraday rotations have been pro- 
duced for the geographic coordinates of the LOFAR core and 
the LOS towards the supernova remnant Cassiopeia A (Cas A; 
RA = 23 h 23 m 27.9 s , DEC = +58°48'42.4"). Cas A was selected 
because it is circumpolar as viewed from LOFAR (its minimum 
elevation is ~ 22°). Therefore, ionFR can calculate the variation 
of (pi on for an entire day. 

Figure [4] shows the variation in the modeled (p mn for two sep- 
arate days: April 1 1th, 2009, close to the most recent minimum 
in solar activity, and April 11th, 2011, by which time solar ac- 
tivity had increase significantly. The next solar maximum is ex- 
pected to be reached around mid 2013. These plots show the 
daily variation of (p lon , increasing around sunrise and decreas- 
ing around sunset. This is expected due to the increase in the 
density of free electrons in the ionosphere from solar irradiation 
during the day. During night hours these free electrons begin to 
recombine with free ions. The degree of variation in 0; on is no- 
ticeably different on the two days. Near solar minimum (Fig. 4] 
left), the daily peak </>j on value is ~ 1 rad irr 2 , while for the 
prediction two years later (Fig. [4] right), in which solar activ- 
ity has increased, the level sometimes exceeds 2 rad irr 2 . Note 
that these time-variable contributions are much larger than the 
formal uncertainties that are achievable through RM-synthesis 
techniques applied to low-frequency radio data (see §5). 

Figure [5] (left) shows weekly averages of the daily maxi- 
mum and minimum modeled |0; on | since April 1998, along the 
LOS of Cas A as seen from the LOFAR core. According to 
the Space Weather Prediction Center (SWPC), Solar Cycle 23 
started in May 1996 and ended in December 2008. ionFR was 
run for each day in Solar Cycles 23 and 24, until April 2012. 
The code produced daily files, each containing 24 values of <pi on . 
From each file the daily maximum and minimum |0j on | values 
were obtained. These values were averaged to give a represen- 
tative maximum and minimum \(/>i n\ for every week since May 
1996. Due to the incompleteness of GIMs within the IONEX 
files for the years 1996, 1997, and part of 1998, the ionospheric 
predictions are only shown since April 1998. The oscillation in 



the minimum |0; on | is a well-known seasonal effect; more ioniza- 
tion is expected in summertime than in winter. In contrast, the 
maximum \<f>i on \ reveals that during the years of greatest solar ac- 
tivity (as reported by SWPC) several rad irr 2 can be reached, as 
viewed from the LOFAR core. It is also noted that when solar 
activity is at its highest, the maximum \<pi on \ no longer appears to 
be dominated by seasonal variations. 

It is evident from Fig. [2] that the ionosphere above the two 
sites chosen for the SKA are subject to the Equatorial Ionization 
Anomaly (EIA). These two regions of enhanced plasma den- 
sity are located approximately 15 degrees north and south of the 
magnetic dip equator ( McDonald et al.|20lT) and are the result 
of the equatorial fountain effect ( Appleton|1946| l. We note that 
the SKA will suffer far higher levels of ionospheric Faraday ro- 
tation than at the LOFAR sites, as the southern component of 
the EIA can pass directly above both locations. Figure [5] (right) 
illustrates this by showing a representative history of |0j on ] for 
the two proposed SKA core sites. Predictions for the two sites 
were made towards the same astronomical object and then av- 
eraged. The LOS chosen to generate the two predictions was 
towards EtaCarinae (RA = 10 h 45 m 03.6 s , DEC = -59°41'04.0"), 
which is circumpolar as viewed from both proposed sites. As 
expected, we observe that |0j on | can be much higher than for 
LOFAR (Fig. 5] left). The seasonal effect is clearly visible in 
the maximum 0; on | curve. Additionally, Figure|2]shows that dur- 
ing 12 hours the TEC and hence ionospheric RM were subject to 
large variations (> 80 TECU). These facts underline the vital im- 
portance of ionospheric calibration for polarimetric studies with 
the SKA — particularly during the day, but even at night. 

Also, as seen in Fig. [2} the EIA passes very close to the site 
of the GMRT, located near Pune, India. Hence, polarimetric ob- 
servations with the GMRT can also benefit greatly from ionFR, 
especially for achieving the full potential of its available bands 
below 300MHz. 

4. RM-synthesis 

RM, which quantifies the amount of Faraday rotation along the 
LOS, has commonly been determined as the gradient of the po- 
larization angle (x) as a function of wavelength squared (A 2 ); 
see, e.g.,|Cooper & Price| ( fl9l^ , |Rand & Lyne| ( fT9S>4"| i, and[Van] 
|Eck et aL ( 201 1] >. These previous RM measure ments assumed 
therefore a linear relationship between ;f and A 2 . Brentjens & de 
|Bruyn| ( |2005| l further developed an alternative method to mea- 
sure Faraday rotation called RM-synthesis, which was first pro- 
posed by Burn (|1966 i and is now increasingly used for such 
measurements (He ald et al.||2009{ |Pizzo et al.||2011| e.g.,). The 
benefits of using RM-synthesis are numerous. For instance, it 
minimizes or eliminates any n7r ambiguity, as opposed to the 
'gradient' method in which x can rotate by 180° an arbitrary 
number of n times between data points at each A 2 . It also uses the 
polarization information across the entire observing bandwidth 
simultaneously such that it is not necessary to detect polarization 
angles at each A 2 . Additionally, it does not assume that there is 
a single RM towards the LOS, i.e. that the linear relationship 
between ;f and A 2 holds. 

Following Brentjens & de Bruyn ( 2005 \ we define the ob- 
served complex polarization vector P (Q + iU) as: 

F((f>)exp[2i(f>A 2 ]d<p, (8) 

oo 

where F (<f>) is the intrinsic complex polarized surface brightness 
per unit Faraday depth, known as the Faraday dispersion func- 
tion (FDF, or Faraday spectrum; |Burn|1966| l. Hence, Eq. [8] can 
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Fig. 4. Prediction of the ionospheric Faraday rotation for two different epochs, as viewed from the LOFAR core along the LOS 
to Cas A. Left: Prediction when solar activity is near its minimum. Right: Prediction when the solar activity rises towards a new 
maximum, expected around May 2013. The dashed vertical lines mark the times of sunrise and sunset. 
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Fig. 5. Weekly averages of the maximum and minimum (blue and red lines, respectively) absolute ionospheric Faraday depth \(pi 0D \ 
from April 1998 - 2012, as modeled by ionFR. The absolute value is shown because the ionospheric Faraday depth is positive or 
negative if observing from the northern or southern hemisphere, respectively. The dashed vertical lines represent the end of Solar 
Cycle 23 (December 2008) and beginning of the current Solar Cycle 24. Left: towards CasA, as viewed from LOFAR. Right: towards 
Eta Carinae, as viewed from an average of the SKA core sites in Western Australia and South Africa. 



be inverted to obtain the FDF from the observed complex polar- 
ization vector: 

P(^ 2 )exp[-2/^ 2 ]d^ 2 . (9) 

oo 

Polarization observations hav e a limited range in A 2 and 
A 2 < is not possible. Therefore, Brentjens & de Bruyn ( 2005 1 



introduce a window function that is non-zero for all observed A 
and zero otherwise. In practice, the integral in Eq.|9]is performed 
as a discrete sum such that for each discrete channel i at A 2 : 



F(4>) 



J^cot ^w,P(^ 2 )exp[-2^(i 2 -^)], (10) 



where is the weight of channel i and A$ is the weighted av- 
erage of the observed A 2 . Considering this inversion as the 'co- 
herent' addition of the observed polarization vectors (a polar- 
ization vector per channel) for a range of Faraday depths, the 
vectors will constructively interfere within the bandwidth for the 
Faraday depth of the observed source and will result in a peak at 
this <p in the FDF. To illustrate this, example FDFs obtained from 
LOFAR High-Band Antenna (HBA) and Low-Band Antenna 
(LBA) observations (see §5) are displayed in Fig. [6] Equation 10 
is used to determine the FDF for the observations described here. 
The relationship between the input observed complex polariza- 
tion vector and the output FDF is given by the rotation measure 
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spread function (RMSF): 



R(<f>) 



Y^oJi ^,exp [-2i<f>(A 2 -A 2 )]. (11) 



RM-synthesis is analogous to performing aperture synthesis 
imaging with an interferometer in the sense that both methods 
make use of Fourier transforms using discrete sampling in A 2 
and spatial frequency coordinates (u, v, w), respectively ( |Heald| 
et al. 12009) e.g.,). Hence, the RMSF in Faraday space </> is anal- 
ogous to t he dirty be am in angular coordinates on the sky {I, m) 
(Bell & EnBlin 20121. As such, fewer gaps in the sampling of A 1 
reduce the side lobes in the RMSF and using larger bandwidths 
in A 2 space increases the resolution in <p space. The resolution 
in Faraday space can be quantified by the FWHM of the RMSF 
function, 

3.S 



FWHM 



rad m 



-2 



(12) 



where A m i n and A max are the shortest and longest observed wave- 
lengths and the constant 3.8 is from Schnitzele r et al.| (|2009). 
Moreover, the uncertainty associated with locating the peak 
Faraday depth in the FDF can be determined in the same way 
as the uncertainty associated with locating the peak flux in an 
aperture synthesis image ( |Fomalont|1999 see): 



FWHM 

2xS/N' 



(13) 



where FWHM is the FWHM of the RMSF defined in Eq.pjand 
S/N is the total polarized signal-to-noise ratio. Equation 13 is 



used to determine the error in the Faraday depth measurements 
presented here. 

5. Model comparison with observational data 

5.1. LOFAR pulsar data 

To measure Faraday depth variations in the ionosphere and to 
compare them with those modeled by ionFR, four bright polar- 
ized pulsars were observed using LOFAR (see Stap pers et al.| 
201 1 1 for a description of LOFAR's pulsar observing modes). 



One or more pulsars were observed on four separate epochs, in- 
cluding times when the ionosphere was expected to be partic- 
ularly dynamic (around sunrise and sunset). See Table [T] for a 
summary of these four observing campaigns. 

The first campaign used the coherently combined LOFAR 
'Superterp j^] to observe PSR B0834+06 in 7 x 10-min integra- 
tions spaced every 50 min. These started 1.8 hr before sunset 
(18:25 UT) and continued until more than 2 hr after astronomi- 
cal twilight (20:40 UT). 

The second Superterp campaign observed PSR B0834+06 in 
20 x 3-min integrations spaced every 7 min. These started after 
astronomical twilight (04:11 UT) and continued until 1.5 hours 
after sunrise (06:08 UT). 

The third Superterp campaign observed PSRs B 1642-03, 
B1919+21 and B2217+47 by cycling consecutively through the 
three pulsars so that each was observed for 12 x 3-min inte- 
grations spaced every 20 min. This enabled quasi-simultaneous 
measurements of the ionospheric Faraday depth variations to- 
wards three widely separated LOSs. These observations started 
before nautical twilight (04:16 UT) and ended two hours after 
sunrise (05:29 UT). 



The fourth campaign observed PSR B0834+06 using the 
LOFAR Superterp stations and two international stations located 
near Nancay, France and near Onsala, Sweden. The pulsar was 
quasi-simultaneously observed by each station in 1 1 x 3-min in- 
tegrations spaced by 17 min. These were done during midday 
when the absolute TEC was expected to be relatively high. This 
enabled measurements of Faraday depth variations from three 
locations separated by long geographical baselines - 594 km 
minimum and 1294 km maximum distance, respectively. 



In all cases, data were written as complex values for the two 
orthogonal linear polarizations. The data were recorded using 
the 200 MHz clock mode, which provides multiple 195.3125 
kHz subbands that are further channelized by an online poly- 
phase filter to 12.2 kHz channels with a time resolution of 81.92 
|xs. Due to limitations on the data rate at the time of observa- 
tion, 6-11 MHz of bandwidth were recorded. In comparison, 
80 MHz of bandwidth can now be recorded by LOFAR in this 
mode. Nonetheless, given the low central observing frequencies 
(~ 125 MHz) the recorded bandwidths were still more than ad- 
equate to achieve precise RM measurements. The complex val- 
ues were converted to 8-bit samples offline and then coherently 
dedispersed and folded using the dspsr program (van Straten & 
Baile s|20TT) l. Radio frequency interference (RFI) was removed 
using the pazi program of PSRCHIVE (|Hotan et al.|2004]>. 



The reduced data were analyzed using an RM-synthesis pro- 
gram in order to determine a precise Faraday depth for each in- 
dividual observation. For each data set, the Stokes parameters 
(IQUV) and associated uncertainties for each frequency chan- 
nel were output for the pulsed section of the pulsar profile using 
the PSRCHIVE program rmfit. The frequency, Q and U in- 
formation were used as the input to the RM-synthesis program, 
which calculated the FDF as a discrete sum for Faraday dep ths 
-50 < < 50 rad irT 2 in steps of 0.001 rad irT 2 using Eq. 10 (see 
Fig.[6]for two examples of these). The peak associate d with the 
instrumental DC signal at ~0 rad m -2 Faraday depth ( Geil et al. 
201 1) see, e.g.,) was subtracted from the FDF before determining 



the peak associated with the Faraday depth towards each pulsar 
LOS. This had no effect on the Faraday peaks of the pulsars since 
the known RM values (ATNF pulsar catalog; Manchester et al. 
2003) ) are > 2xFWHM of the RMSF, see Eq. |12| The Faraday 
depth at which the peak in the FDF occurred was assumed to be 
the measured Faraday depth of the ISM and ionosphere towards 
the pulsar, </>ism + 0ic»LJ For each observing campaign, we es- 
timated the instrumental error by measuring the Faraday depth 
of the instrumental polarization peak around rad m~ 2 in the 
FDF of each observation, weighted by the S/N, and taking the 

2|. This demonstrated that larger 
ligher S/Ns also tended to reduce 
the scatter in instrumental Faraday depth around rad m -2 . The 
total error on the Faraday depth was taken to be the formal error 
from Eq. 13 added in quadrature with the instrumental error. The 
linear polarization had S/N > 30 for all observations; this is well 
above the threshold necessary for reliable Faraday depth mea- 
surements (see Macquart et al. 2012 and references therein). 



rad m 2 +lcr limits (Table 
bandwidth observations with 



While we have shown that LOFAR provides reliable Faraday 
depths, we note that absolute polarization calibration (e.g. to de- 
termine absolute polarization angles) has not yet been applied to 
the data. 



10 This assumes no Faraday rotation in the pulsar magnetosphere it- 
self; see Noutsos et al. 2009 and Wang et al.|2011| for a discussion on 



The 330-meter-wide inner core of the array, which hosts 6 stations. possible Faraday rotation within pulsar magnetospheres. 
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Fig. 6. Examples of the absolute value of the normalized Faraday Dispersion Function (FDF), obtained from two LOFAR observa- 
tions. The peaks at zero Faraday depth, F{<p ~ 0), indicate the instrumental response. Left: FDF obtained from a single 3-minute 
HBA observation (119 - 125 MHz) of PSR B2217+47. The peak at F(4> = -34.08 rad irT 2 ) is the response due to the polarized 
flux of the pulsar. The FWHM of the RMSF is 6.6 rad rrr 2 . Right: FDF obtained from a single 3-minute LBA observation (58 - 64 
MHz) of PSR B1919+21. The polarized flux of the pulsar is responsible for the peak at F(<p = -15.94 rad m~ 2 ). The FWHM of the 
RMSF is 0.84 rad m- 2 . 



No. 


PSR 


Date 


Obs. duration 


Surprise, set] time 


LOFAR stations 


Time obs. 


No. obs. 


Freq 


Elevation 


LOFAR obs. IDs 




(B-name) 


(dd.mm.yyyy) 


(hh:mm UT) 


(hh:mm UT) 




(min) 




(MHz) 


(deg) 




1 


B0834+06 


11.04.2011 


16:40-22:50 


18:25 


CS00[2-7]HBA 


10 


7 


120-126 


20-45 


L25152-L25158 


2 


B0834+06 


20.10.2011 


04:20-07:33 


06:08 


CS00[2-7]HBA 


3 


20 


129-140 


37-44 


L32350-L32369 


3 


B1642-03 


23.03.2012 


04:08 - 07:40 


05:29 


CS00[2-7]HBA 


3 


12 


119-125 


16-34 


L53966-L53977 


3 


B1919+21 


23.03.2012 


04:04-07:36 


05:29 


CS00[2-7]LBA 


3 


12 


58-64 


45-59 


L53942-L53953 


3 


B2217+47 


23.03.2012 


04:12-07:44 


05:29 


CS00[2-7]HBA 


3 


12 


119-125 


37-72 


L53990-L54001 


4 


B0834+06 


10.07.2012 


11:20-14:43 


N/A 


CS00[2,3,5-7]HBA 


3 


1 1 


119-129 


38-44 


L6 1 473 -L6 1483 


4 


B0834+06 


10.07.2012 


11:25-14:48 


N/A 


FR606HBA 


3 


1 1 


119-129 


42-49 


L61532-L61542 


4 


B0834+06 


10.07.2012 


11:30-14:53 


N/A 


SE607HBA 


3 


1 1 


119-129 


32 - 39 


L61520-L61530 



Table 1. Summary of the four LOFAR pulsar observing campaigns. Columns 1-1 1 show the number of the observing campaign, 
pulsar B names, date, duration, time of sunset or sunrise during the observations, LOFAR station(s) (CS00[2-7] indicate the LOFAR 
'Superterp' stations combined in tied-array mode, HBA and LBA indicate High-Band Antenna stations and Low-Band Antenna 
stations, respectively), individual observation integration times, total number of observations, frequency range, elevation range, and 
LOFAR observation identification numbers (obs. IDs). 



The observed and ionFR-modeled Faraday depths as a func- 
tion of time for the four LOFAR observing campaigns are plot- 
ted in Figs. [7] [8] and [9] In general, the modeled and observed 
Faraday depth variations agree very well. However, there are a 
few instances where the observed and modeled values differ by 
more than lcr. The measured RMs could still be affected by in- 
terference in some cases and it is also quite possible that there are 
unmodeled ionospheric variations on short timescales (see also 
Sj6}. TEC data from CODE was used for the B0834+06 sunset 
and sunrise campaigns. CODE and ROB TEC data were both 
available for the observations which took place in 2012 and are 
compared in Figs. [8] and [9] In each case, there is a clear trend 
in the Faraday depth as a function of time. At sunrise and sun- 
set there is a particularly distinct variation in Faraday depth of 
approximately 2 rad nT 2 , as expected (cf. Fig. [4]). Even during 
midday, Fig. [9] when the absolute TEC is expected to be high 
but relatively constant (e.g.,|Kassim et al. 2007 Fig. 14), smaller 



variations in Faraday depth are still evident and well modeled in 
general. 



Figs. [8] and [9] compare the modeled Faraday depths using 
CODE and ROB TEC data. In both cases, the observations and 
model show good agreement, although it is evident that there 
are variations on timescales less than one hour which are not 
resolved in the model using CODE data. The model using ROB 
TEC data more accurately fits the observations. The finer griding 
and smaller RMS available with the ROB TEC data also allows 
smaller fluctuations in Faraday depth to be resolved. For the ob- 
servations of B0834+06 using the Superterp and international 
stations over midday, Pig. [9j this is especially significant, where 
the variations in the data appear on shorter timescales and in 
Faraday depths within the error bars of the model using CODE 
data. The model output of ionFR shows good agreement with 
the measured Faraday depths for all of the LOFAR observing 
campaigns depicted in Figs. [7] [8] and [9] Therefore, this model 
was used to subtract the contribution of the ionospheric Faraday 
depth from the measurements in order to determine the Faraday 
depth of the ISM, 0ism, in the direction of these pulsars. 

To determine 0ism for each pulsar, the constant offset which 
yielded the minimum weighted chi-squared value between the 
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Fig. 7. Observed Faraday depths, (^observed, and ionFR-modeled ionospheric Faraday depths, (p lon , towards PSR B0834+06 as a 
function of time (observations: blue circles, left axis labels; model: red triangles, right axis labels). Only CODE TEC data was 
available for these predictions. Left panel: seven LOFAR Superterp HBA observations during sunset. Right panel: twenty LOFAR 
Superterp HBA observations during sunrise. The vertical dashed lines indicate the time of sunset and sunrise, respectively. The 
offsets between ^observed an d 0k>n for the sunset and sunrise campaigns are 25.15 ±0.18 and 24.94 ± 0.24 rad irT 2 , respectively. 



pulsar Faraday depth measurements and the ionospheric mod- 
eled Faraday depths from ionFR was used (see Table [2}. Using 
the CODE data, the reduced chi-squared values range from 0.4 
- 1 .0, whereas with the ROB data the reduced chi-squared val- 
ues range from 0.9 - 1.3. The residual differences between the 
observations and model may be due to small scale variations 
in the ionosphere which affect the observations but which are 
not resolved due to the time resolution of the TEC data. The 
standard deviations of (^observed - <^ion divided by the square root 
of the number of measurements in each observational campaign 
range from 0.03 - 0.09 rad irr 2 for the CODE output and from 
0.02 - 0.07 rad rrT 2 for the ROB output. This indicates that 
both fit the measurements well, although the ROB data gives 
reduced chi-squared values closer to 1 and smaller errors due to 
the smaller uncertainties compared with the CODE maps. Also, 
both sources of TEC maps give consistent results, where all 0ism 
values obtained using CODE and ROB data agree within l-cr for 
the same observing campaign and LOS. The Faraday depth val- 
ues for B0834+06 obtained from the three observational cam- 
paigns over six months apart also agree at or below the 2-cr level 
for both CODE and ROB data. The weighted mean for all CODE 
and ROB data available for B0834+06 gives ISM = 25.12 + 0.07 
rad m" 2 (5 LOSs) and ISM = 25.26 + 0.05 rad irr 2 (3 LOSs), 
respectively. Although there are two more data sets with CODE 
data available, the errors assumed for the ROB data are smaller. 
It is worth noting that these are among the most precise pulsar- 
derived 4>\sm measurements ever obtained ( [Manchester et al7| 
|2005[ see ATNF catalog), and improve significantly on the pre- 
cision of previous RM measurements for these pulsars. 

In order to determine the Faraday depth of B0834+06 us- 
ing an independent instrument, on April 22nd, 2012 at 16:50 UT 
we performed a 20-minute observation of PSR B0834+06 using 



the WSRT with the PuMall pulsar backend (Karuppusamy et al. 
2008 ) from 3 10 to 390 MHz. We recorded baseband data, which 
were folded, dedispersed and subsequently analyzed using the 
same RM-synthesis method described above, yielding ^observed = 
28.2 ± 1.8 rad nr 2 . The ionFR code predicts (f> ion = 2.75 + 0.15 
rad irr 2 , using CODE maps, for the given time and LOS of this 



observation, resulting in a corrected value of (f>^f T = 25.4 ± 1.8 
rad irr 2 . Using ROB maps we calculated a corrected value of 
^iSM* 7 = 25.1 + 1.5 rad irr 2 . These values are in excellent agree- 
ment with the more precise value derived from the LOFAR ob- 
servations. The precision of the WSRT measurement is lower in 
part because of the higher observing frequency, see Eqs.[T2]and 



13 This demonstrates the power of low-frequency observations 
for the purpose of determining accurate 0ism- 

Our measured Faraday depths are consistent with previously 
published measurements for these pulsars, but are significantly 
more precise. The values of 0ism obtained from the LOFAR ob- 
servations of B0834+06 agree to within 2.5 cr of the catalog 
value for this pulsar, RM psrcat = 23.6 ± 0.7 rad irr 2 . It is un- 
clear, but likely, that the catalog value was calibrated for iono- 
spheric Faraday rotation (see Hamilton & Lyne 1987, §2). That 
paper also states that the ionosphere contributes 1-8 rad m~ 2 
to other observations, using either a geostationary satellite or 
an ionosonde, and that the subtraction of this introduces un- 
certainties of approximately 1-2 rad irr 2 . The 2.5-cr difference 
between the 0ism obtained using LOFAR and the ATNF cata- 
log value determined in 1987, plus the possibility that pulsar 
RMs may change on multi-year timescales ( Weisberg et al. 2004 
e.g.,) due to variations in the polarized pulsar emission and/or 
electron density changes along the LOS through the ISM, pro- 
vided the motivation for the recent independent comparison ob- 
servation of PSR B0834+06 using the WSRT. 

The values of </>ism obtained for the three pulsars observed 
quasi-simultaneously are also very precise and are in excellent 
agreement with those of the ATNF pulsar catalog. Prior mea- 
surements for B1642-03 and B1919+21, also calibrated for 
ionospheric Faraday rotation using geostationary satellite and 
ionosonde data, give RMism = 15. 8 + 0.3 rad irr 2 and RMj sm = 
-16.5 + 0.5 rad m" 2 , respectively (Ha milton & Lyne|l987 ), and 
are also in good agreement within 1 cr of the values obtained 
in this work (ROB-subtracted values 4>\su = 16.04 + 0.18 rad 
irr 2 and 4>i$ M = -16.92 + 0.07 rad rrT 2 ). A prior measurement 
for B2217+47, RM ISM = -35.3 ± 1.8 rad irr 2 , also calibrated 
for ionospheric Faraday rotation using a geostationary satellite 
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Fig. 8. Observed Faraday depths, o b se rved> an d ionFR-modeled ionospheric Faraday depths, 0; on , towards three pulsars as a function 
of time during sunrise (observations: blue circles, left axis labels; model: red triangles, right axis labels). Upper panels: twelve 
LOFAR Superterp HBA observations of PSR B 1642-03. Middle panels: twelve LOFAR Superterp LBA observations of PSR 
B1919+21. Lower panels: twelve LOFAR Superterp HBA observations of PSR B2217+47. The vertical dashed lines indicate the 
time of sunrise, (a) shows the ionFR model using CODE TEC data and IGRF11. (b) shows the ionFR model using ROB TEC data 
andlGRFll. 



( |Manchester|1972| ), is also in good agreement within 1 cr of the 
value derived here (ROB-subtracted value 0ism = -35.60 ± 0.1 1 
rad irT 2 ). 

Together these observations provide a convincing verifica- 
tion of the accuracy of the ionFR model and demonstrate the 
ability to derive Faraday depths resulting solely from the ISM 
by robustly removing the time and direction-dependent Faraday 
depth introduced by the ionosphere, even during some of the 
most turbulent periods expected in daily ionospheric variations; 
see §6 for further discussion. 



for each slot. From the polarization maps we then calculated the 
polarization angle of PSR B1937+21. 



Fig. 



10 also shows the predicted rotation due to the iono- 

using CODE 
and then 



sphere. This was calculated by modeling the 
TEC data (no ROB data is available prior to 2012) 
converting the ionospheric Faraday depths into polarization an- 
gles at 118.2 MHz. After removing a constant offset, the mod- 
eled Faraday depths and observation match very well. This com- 
parison provides further verification of ionFR-modeled Faraday 
depths and also shows that the ionospheric correction can be 
done for interferometric imaging data. 



5.2. WSRT imaging data 

Here we compare the ionFR-modeled Faraday depths with 
archival 1 18.2-MHz WSRT observations of PSR B 1937+21. 

PSR B1937+21 was observed with the WSRT Low 
Frequency Front Ends on 25 March 2005 for 3 hours following 
sunrise to monitor ionospheric Faraday rotation. The observed 



polarization angle at 118.2 MHz is shown in Fig. 10 and it is 



seen to increase rapidly and smoothly. The polarization angle 
presented refers only to a 4 MHz band, from 116.05 - 120.34 
MHz. The calibrated data were split into 30 time slots, of 6 min- 
utes each, in order to produce a pair of Stokes Q and U maps 



As this example shows, low-frequency observations lasting a 
few tens of minutes or more suffer depolarization from variable 
ionospheric Faraday rotation; in this case, the polarization angle 
has made almost a full turn within 3hrs. The time-variable iono- 
spheric Faraday depth will lead to misalignment between the po- 
larization vectors, which therefore results in significant signal 
loss. Even an instantaneous measurement gives a Faraday depth 
that is systematically biased by ±0.25 - 2.5 rad rrT 2 or more, 
which can be a large fraction of the total measured Faraday depth 
for nearby sources (d < 2 kpc), greatly impeding our ability to 
use these as probes of the local interstellar magnetic field. 
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Fig. 9. Observed Faraday depths, ^observed) an d ionFR-modeled ionospheric Faraday depths, <p mn , towards B0834+06 as a function 
of time during midday (observations: blue circles, left axis labels; model: red triangles, right axis labels). Upper panels: eleven 
LOFAR HBA observations using the international station near Onsala, Sweden. Middle panels: eleven LOFAR Superterp HBA 
observations. Lower panels: eleven LOFAR HBA observations using the international station near Nancay, France, (a) shows the 
ionFR model using CODE TEC data and IGRF11. (b) shows the ionFR model using ROB TEC data and IGRF11. 



25/03/2005: CODE 



I "S I-: I ■ ■: .7 21 

TEc i iGRFn 



: |Homtm| 



Fig. 10. PSR B 1937 +21 's observed polarization angle compared 
with the predicted amount of rotation due to the ionosphere (red 
and blue points, respectively). The errors on the polarization an- 
gles for the WSRT data are less than 2 degrees and are not visible 
in the plot. The vertical dashed line marks the time of sunrise. 



6. Discussion 

The Faraday depths of the ISM towards PSRs B0834+06, 
B1642-03, B1919+21 and B2217+47 presented here (Tableg 
are some of the first determined using the RM-synthesis tech- 
nique and are among the most precise Faraday depth measure- 
ments ever made towards a pulsar. In the ATNF catalog, there 
are only a dozen measurements with absolute precision < 0.2 
rad irT 2 , determined for some of the brightest known pulsars. 
The fact that the Faraday depths obtained here are also in excel- 
lent agreement with the catalog and/or alternatively determined 
values also demonstrates the robustness of ionFR for subtract- 
ing the 0; on contribution and hence for allowing us to reap the 
benefits of low-frequency measurements. 

Using multi -beaming capabilities (Stapp ers et al.||2011| l, 
LOFAR will provide high-precision Faraday depth measure- 
ments for hundreds of nearby pulsars, which can be used to 
probe the interstellar magnetic field of the Galaxy (e.g. Sobey 
et al. 2013, in preparation). Since the observations presented 
here used just a tenth of the now-available LOFAR bandwidth, 
the prospects for even higher-precision Faraday depth measure- 
ments are excellent, especially using LOFAR LBA observations 
(10 - 90 MHz). However, with the TEC data available from 
ROB these will be limited to precisions of approximately 0.05 
rad irT 2 because of the systematic uncertainty provided by the 
RMS VTEC maps. In other words, more sophisticated calibra- 
tion techniques will need to be devised in order to reap the 
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PSR 


Date 


Telescope 


Station(s) 


FWHMhsf 




.CODE 


y 2 

™ red 


.ROB 
t^ISM 


y 2 

A red 


ISM 




[dd.mm.yy] 






[rad nT 2 ] 


[rad nT 2 ] 


[rad m" 2 ] 


CODE 


[rad nT 2 ] 


ROB 


[rad m ] 


B0834+06 


1 1.04.1 1 


LOFAR 


Superterp 


6.2 


0.20 


25.15 ±0.18 


0.52 


N/A 


N/A 


23.6 ± 0.7 


B0834+06 


20.10.11 


LOFAR 


Superterp 


9.4 


0.20 


24.94 ± 0.24 


0.64 


N/A 


N/A 


23.6 ± 0.7 


B 1642-03 


23.03.12 


LOFAR 


Superterp 


6.6 


0.20 


15.98 ± 0.23 


0.8 


16.04 ± 0.18 


1.3 


15.8 ± 0.3 


1_> 1717Til 


23.03.12 


LOFAR 


Siinprtpm 


0.8 


0.08 


-16 qs + n i? 

1U.7J IT_ U.lz. 


0.5 


- 1 6 9? + 07 


1.3 


-16.5 + 0.5 


B22 17+47 


23.03.12 


LOFAR 


Superterp 


6.6 


0.13 


-35.72 + 0.15 


1.0 


-35.60+ 0.11 


1.1 


-35.3 ± 1.8 


B0834+06 


22.04.12 


WSRT 


N/A 


12.0 


1.2 


25.4 ± 1.8 


N/A 


25.1 ± 1.5 


N/A 


23.6 ± 0.7 


B0834+06 


10.07.12 


LOFAR 


Superterp 


4.2 


0.07 


25.16 + 0.13 


0.5 


25.23 ± 0.08 


0.9 


23.6 ± 0.7 


B0834+06 


10.07.12 


LOFAR 


FR606 


4.2 


0.10 


25.21 ± 0.15 


0.4 


25.16 + 0.10 


1.2 


23.6 ± 0.7 


B0834+06 


10.07.12 


LOFAR 


SE607 


4.2 


0.15 


25.22 + 0.18 


0.6 


25.39 + 0.14 


1.0 


23.6 ± 0.7 



Table 2. Summary of results from the four LOFAR pulsar observing campaigns and WSRT observation. Columns 1-1 1 show the 
pulsar observed, date (see Table [TJ for specific times), the telescope, the LOFAR stations (if applicable, see Table [TJ for specific 
station names), the FWHM of the RMSF from RM-synthesis, the error introduced by instrumental effects <x inst , the Faraday depth 
of the ISM towards the pulsar ^y&m as determined using the CODE TEC maps, the reduced chi-squared value x] eA f° r comparison 
between the ionospheric Faraday depths from ionFR using CODE TEC data and the observed Faraday depths, the Faraday depth of 
the ISM towards the pulsar </>ism as determined using the ROB TEC maps, the reduced chi-squared value for comparison between 
the ionospheric Faraday depth produced from ionFR using CODE TEC data and the observed Faraday depths. The final column 
gives the catalog RM value from psrcat (Manchester et al. 2005| for comparison. Note that ROB TEC maps were not available for 
the observations of B0834+06 taken in 201 1. Since a single observation was obtained with WSRT, the reduced chi-squared value is 
not applicable because the Faraday depth of the ISM was determined by subtracting the Faraday depth of the ionosphere from the 
observed Faraday depth. 



full benefit of LOFAR's low observing frequencies and large 
fractional bandwidth. For instance, it is possible to measure 
time-dependent, differential TEC and Faraday rotation between 
LOFAR stations using visibility data from observations of bright 
(even unpolarised) calibrators. These data can be combined with 
magnetic field models in order to predict the absolute Faraday 
depth, in addition to using the TEC data described here for com- 
parison. This is currently being investigated using LOFAR long- 
baseline observations, but will be less suitable for more compact 
array designs (e.g., LOFAR pulsar observations typically use 
only the 2-km core of the array). Higher-precision Faraday depth 
measurements are especially important for determining possible 
long-term variations of the Faraday depth of pulsars due to fluc- 



tuations in their magnetosphere or the ISM (|Weisberg et al. 2004 



e.g.,), particularly as some pulsars have large relative velocities. 

Calibrating for 0; on is also important for higher-frequency 
observations. Assuming a bandwidth of 1.21 - 1.51 GHz, such 
as that of the multibeam receiver used at the 100-m Effelsberg 
radio telescope (Ba rr et al.|2013[ ) , the theoretical FWHM of the 
RMSF is 173 rad nT" 1 . This yields an uncertainty in Faraday 
depth of 8.7 rad nT 2 given a S/N of 10, which is comparable to 
the 0j on reached during solar maximum as Fig. [5] (left) shows. 
Given larger bandwidths, such as those planned for the SKA and 
its pathfinders, this problem becomes worse and is relevant even 
during night-time and solar minimum. It is also clear that correc- 
tions for the ionospheric Faraday depth are important for obser- 
vations of pulsars and extragalactic sources, particularly towards 
the halo of the Galaxy, since the Faraday rotation expected is of- 
ten much lower than those located towards the plane. For the 650 
pulsars with rotation measure data, 284 of these are located over 
300 pc above or below the galactic plane and have a median ro- 
tation measure of ~39 rad nT 2 . Ionospheric Faraday depth can 
thus contribute significantly to the total observed Faraday depth, 
especially at times near solar maximum or when observing from 
lower latitude sites. 

The TEC maps from ROB will help to perform differen- 
tial Faraday rotation studies between LOFAR core stations and 
the international stations. Fig.|3]shows a VTEC map from ROB 
over Europe. It is clear that the international stations are subject 



to different ionospheric conditions than the core, which results 
in different amounts of Faraday rotation for the signals arriv- 
ing at each station. This effect needs to be calibrated for be- 
fore combining all the LOFAR stations to carry out polariza- 
tion studies. On the other hand, to better understand the iono- 
spheric variations above just the Dutch LOFAR stations (max. 
baseline < 100 km), better geographically resolved TEC maps 
are needed. Alternatively, raw, dual-frequency GPS data should 
be directly analyzed. Raw GPS data is also desirable to use 
when short timescale ionospheric changes (on the order of sec- 
onds to a few minutes) need to be studied. For arrays located in 
Europe this is possible due to the considerable number of GPS 
stations from the permanent European network. However, this 
is not yet the case for arrays like the GMRT or future SKA. 
For instance, in India there is only one active station [^] (near 
Bangalore) which is part of the IGS network. This is also the 
case for South Africa and Western Australia, where there are a 
few more stations but they are very far apart. Regional dedicated 
networks, like in Europe, are needed for these arrays in order 
to gain a more detailed picture of the ionospheric sky over these 
sites. Alternatively, high accurate ionospheric corrections for the 
SKA could be possible through its own imaging data; a route be- 
ing explored, as previously mentioned, for LOFAR. 

A single LBA station operating at 20 MHz has a beam width 
of ~ 13°, which spans less than 2 cell widths in the ROB TEC 
maps. This is because each of the cells (0.5°) extends ~ 60 km at 
the altitude of the ionospheric thin shell, which corresponds to 
a resolution of ~ 7.5° in the plane of the sky. When using ROB 
data, this implies that we expect to predict the same ionospheric 
Faraday rotation for any source located within a radius of 7.5°. 

We are expanding ionFR to support also the use of TEC data 
from the last release of the IRI (IRI-2007^- This addition will 
allow access to TEC data from before 1995. This type of data 
was not included in the analysis herein due to the fact that no 
uncertainties are available for the IRI. 



See http : //ig scb . j pi . nasa . gov/ network/hourly . html 
12 See http : // omniweb . gsf c . nasa . gov/vitmo/ir i_vitmo . html 
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PSR B0834+ 06 (CODE) 
PSR B08S4 + 06 (ROB) 



Si 1 1 kt1 i T] ■ 



I ll III I" X 

Fig. 11. The measured Faraday depth of the ISM, 0ism> towards 
PSR B0834+06, as determined from HBA observations using 
the LOFAR FR606, Superterp and SE607 stations. These are 
plotted as a function of the geographic latitudes of the stations. 
The data points are offset slightly in latitude for ease of compar- 
ison. The measurements are also compared with the range in air- 
mass calculated for the station location towards PSR B0834+06 
for the duration of the observations (red triangles, right axis). 
The ionFR model was run using both TEC data from CODE 
(green squares) and ROB (blue points). 



To further improve the accuracy of the ionFR- modeled 
Faraday depths, the Earth will be considered as an ellipsoid 
and also the various ionospheric layers (D, E, and F) can be 
treated. A three-dimensional model for the ionosphere should 
also be possible using ray tracing for three-dimensional tomog- 
raphy of the ionosphere (Bust & Mitchell 2008, e.g.,). An indi- 
cation that this may be needed for future high-accuracy Faraday 
depth determinations with LOFAR is the increase in the derived 
Faraday depth of the ISM towards B0834+06 as a function of 
latitude in the international station observations, see Fig. 11 
Since the elevation of the pulsar is lower at higher geographic 
latitudes in these quasi-simultaneous observations, the LOS to- 
wards B0834+06 passes through a larger ionospheric depth as 
seen from Sweden compared with the Netherlands and France. 
Fig. 1 1 shows how the 0ism determined at three latitudes appears 
to follow a similar trend to the airmass. Therefore, a correction 
similar to, although not quite as large as, the airmass may be 
needed (Wielebinski & Shakeshaf t|1962 e.g.,). 



7. Conclusions 

We have presented ionFR, a code that models the ionospheric 
Faraday rotation using publicly available TEC maps and the 
IGRF11. In §3 we show modeled ionospheric Faraday depths 
for changing levels of solar activity and different geographic lo- 
cations. In §5 we compare ionFR-modeled Faraday depths with 
low-frequency data from LOFAR and the WSRT These obser- 
vational comparisons demonstrate the robustness and accuracy 
of the modeled data. 

We have shown in §5.1 that calibrating LOFAR data with 
ionFR provides very high-precision pulsar RMs (absolute er- 
ror < 0.1 rad m -2 ). The applicability of precise RMs is broad. 
For instance, precise RMs can be used to map the structure of 
the Galactic magnetic field. Also, it is now possible to monitor 



pulsar RMs on multi-year timescales. For example, B0834+06, 
shows an apparent increase in its RM after ~ 25 years, assuming 
previous ionospheric calibration was correct. 

Our code represents an alternative, and moreover, a cheaper 
solution when no GPS receivers are co-located with the radio 
telescope carrying out the observation. Additionally, GPS re- 
ceivers may need periodic maintenance, which requires an in- 
vestment of time and money. Therefore, this code is put forward 
as a simple and costless method to the community to predict and 
correct for ionospheric Faraday rotation. 

ionFR will be used to correct for the ionospheric Faraday ro- 
tation in projects such as the WSRT Continuum Legacy Survey 
at 350 MHz: Low-frequency Galaxy Continuum Survey (targets 
were observed in 2012), the LOFAR Magnetism Key Science 
Project ( |Anderson et al.||2012} MKSP;) and the Polarisation 
Sky Survey of the Universe's Magnetism (Gaen sler et al.|2 010 



POSSUM;) project planned with the Australian Kilometre Array 
Pathfinder (Johnsto net al.|2007| ASKAP;) telescope. 

Lastly, Phase I of the SKA will provide an unprecedented 
low-frequency radio telescope, capable of making, e.g., a de- 
tailed map of the Galactic magnetic field structure through 
Faraday depth measurements of pulsars. However, given that the 
ionospheric equatorial anomaly sometimes passes directly over 
the two sites in South Africa and Western Australia, it will be 
crucial to have a robust and accurate calibration procedure in 
place to take full advantage of what the SKA has to offer. 

Acknowledgements. The authors want to express their gratitude to Dr. Dominic 
Schnitzeler and Dr. David Champion for their valuable comments and sugges- 
tions to improve the quality of the paper. CSB and CS would like to thank the 
Deutsche Forschungsgemeinschaft (DFG) which funded this work within the re- 
search unit FOR 1254 "Magnetisation of Interstellar and Intergalactic Media: 
The Prospects of Low-Frequency Radio Observations". JWTH was funded for 
this work by a Veni Fellowship of the Netherlands Foundation for Scientific 
Research (NWO). OW is supported by the DFG (Emmy-Noether Grant WU 
588/1-1) and by the European Commission (European Reintegration Grant 
PERG02-GA-2007-224897 WIDEMAP). CF and GM acknowledge financial 
support by the Agence Nationale de la Recherche through grant ANR-09-JCJC- 
0001-01. The Low Frequency Array (LOFAR) was designed and constructed 
by ASTRON, the Netherlands Institute for Radio Astronomy, and has facili- 
ties in several countries, that are owned by various parties (each with their own 
funding sources), and that are collectively operated by the International LOFAR 
Telescope (ILT) foundation under a joint scientific policy. The WSRT is operated 
by ASTRON/NWO. 



References 

Afraimovich, E. L., Smolkov. G. Y., Tatarinov, P. V., & Yasukevich, Y. V. 2008, 
in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 6936, Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series 

Anderson, J., Beck, R., Bell, M., et al. 2012, ArXiv e-prints 

Appleton, E. V. 1946, Nature, 157, 691 

Barr, E. D., Guillemot, L., Champion, D. I., et al. 2013, MNRAS, 429, 1633 
Beck, R. 2009, in Revista Mexicana de Astronomia y Astrofisica Conference 

Series, Vol. 36, Revista Mexicana de Astronomia y Astrofisica Conference 

Series, 1-8 

Bell, M. R. & Enfilin, T. A. 2012, A&A, 540, A80 

Bilitza, D. & Reinisch, B. W. 2008, Advances in Space Research, 42, 599 

Brentjens, M. A. 2008, A&A, 489, 69 

Brentjens, M. A. & de Bruyn, A. G. 2005, A&A, 441, 1217 

Brown, J. C, Haverkorn, M., Gaensler, B. M., et al. 2007, Apl, 663, 258 

Burn, B. J. 1966, MNRAS, 133, 67 

Bust, G. S. & Mitchell, C. N. 2008, Reviews of Geophysics, 46 
Cooper, B. F. C. & Price, R. M. 1962, Nature, 195, 1084 

Erickson, W. C, Perley, R. A., Flatters, C, & Kassim, N. E. 2001, A&A, 366, 
1071 

Finlay, C. C, Maus, S., Beggan, C. D., et al. 2010, Geophysical Journal 

International, 183, 1216 
Fomalont, E. B. 1999, in Astronomical Society of the Pacific Conference Series, 

Vol. 180, Synthesis Imaging in Radio Astronomy II, ed. G. B. Taylor, 

C. L. Carilli, & R. A. Perley, 301 



13 



C. Sotomayor-Beltran et al: Ionospheric Faraday rotation 

Gaensler, B. M., Landecker, T. L., Taylor, A. R., & POSSUM Collaboration 

2010, Bulletin of the American Astronomical Society, 42, #470.13 
Garrett, M. A., Cordes, J. M., Deboer, D. R., et al. 2010, in ISKAF2010 Science 

Meeting 

Geil, P. M., Gaensler, B. M., & Wyithe, J. S. B. 2011, MNRAS, 418, 516 
Greisen, E. W. 2003, Information Handling in Astronomy - Historical Vistas, 
285, 109 

Hamilton, P. A. & Lyne, A. G. 1987, MNRAS, 224, 1073 
Han, I. L., Manchester, R. N., Lyne, A. G, Qiao, G. J., & van Straten, W. 2006, 
Apl, 642, 868 

Heald, G, Braun, R., & Edmonds, R. 2009, A&A, 503, 409 
Hotan, A. W., van Straten, W., & Manchester, R. N. 2004, PASA, 21, 302 
lohnston, S., Bailes, M., Bartel, N, et al. 2007, PASA, 24, 174 
Karuppusamy, R., Stappers, B., & van Straten, W. 2008, PASP, 120, 191 
Kassim, N., White, S., Rodriquez, P., et al. 2010, in Advanced Maui Optical and 

Space Surveillance Technologies Conference 
Kassim, N. E., Lazio, T. J. W., Erickson, W. C, et al. 2007, ApJS, 172, 686 
Macquart, I. -P., Ekers, R. D., Feain, I., & lohnston-Hollitt, M. 2012, ArXiv e- 

prints 

Manchester, R. N. 1972, Apl, 172, 43 

Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, VizieR Online 

Data Catalog, 7245, 
McDonald, S. E., Coker, C, Dymond, K. E, Anderson, D. N., & Araujo-Pradere, 

E. A. 201 1, Radio Science, 46, RS6004 
Mitchell, D., Greenhill, L. I., Clark, M., et al. 2010, in RFI Mitigation Workshop 
Nigl, A., Zarka, P., Kuijpers, J., et al. 2007, A&A, 471, 1099 
Noutsos, A., lohnston, S., Kramer, M., & Karastergiou, A. 2008, MNRAS, 386, 

1881 

Noutsos, A., Karastergiou, A., Kramer, M., Johnston, S., & Stappers, B. W. 2009, 

MNRAS, 396, 1559 
Pang, B., Pen, U.-L., Matzner, C. D., Green, S. R., & Liebendorfer, M. 2011, 

MNRAS, 415, 1228 
Pizzo, R. F, de Bruyn, A. G, Bernardi, G, & Brentjens, M. A. 2011, A&A, 525, 

A104 

Rand, R. J. & Lyne, A. G. 1994, MNRAS, 268, 497 

Schaer, S. W., Gurtner, W., & Feltens, J. 1998, in Proceedings of the 1998 IGS 

Analysis Centres Workshop, ESOC, Darmstadt, Germany, 233-247 
Schnitzeler, D. H. F. M., Katgert, P., & de Bruyn, A. G. 2009, A&A, 494, 61 1 
Sokoloff, D. D., Bykov, A. A., Shukurov, A., et al. 1998, MNRAS, 299, 189 
Stappers, B. W, Hessels, J. W. T, Alexov, A., et al. 2011, A&A, 530, A80 
Swamp, G. 1991, in Astronomical Society of the Pacific Conference Series, 

Vol. 19, IAU Colloq. 131: Radio Interferometry. Theory, Techniques, and 

Applications, ed. T. J. Comwell & R. A. Perley, 376-380 
Van Eck, C. L., Brown, J. C, Stil, J. M., et al. 201 1, ApJ, 728, 97 
van Straten, W. & Bailes, M. 2011, PASA, 28, 1 
Wang, C, Han, J. L., & Lai, D. 2011, MNRAS, 417, 1183 
Weisberg, J. M., Cordes, J. M., Kuan, B., et al. 2004, ApJS, 150, 317 
Wielebinski, R. & Shakeshaft, J. R. 1962, Nature, 195, 982 



14 



